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Abstract 

The spectral density for two dimensional continuum QCD has a non-analytic be- 
havior for a critical area. Apparently this is not reflected in the Wilson loops. However, 
we show that the existence of a critical area is encoded in the winding Wilson loops: 
Although there is no non-analyticity or phase transition in these Wilson loops, the 
dynamics of these loops consists of two smoothly connected domains separated by the 
critical area, one domain with a confining behavior for large winding Wilson loops, and 
one (below the critical size) where the string tension disappears. We show that this 
can be interpreted in terms of a simple tunneling process between an ordered and a 
disordered state. In view of recent results by Narayanan and Neuberger this tunneling 
may also be relevant for four dimensional QCD. 
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1 Introduction 



Recently there appeared an interesting paper by Narayanan and Neuberger [1] where four 
dimensional Wilson operators for QCD on a lattice were shown to have a finite and nontrivial 
continuum limit. One of the surprising results is that the spectral density of eigenvalues of 
the Wilson loop develops a gap at infinite N for small size loops. This gap is closed when 
the loop size is increased sufficiently The resulting eigenvalue distribution was compared to 
the corresponding density in two dimensional continuum QCD, and the agreement between 
four and two dimensions was very good. 

The spectral density was computed by Durhuus and the author [2] from the Makeenko- 
Migdal equation. Later on the same result was obtained by Bassetto, Griguolo, and Vian 
[3] by direct summation of winding Wilson loops. One result is that for a critical size of 
the (dimensionless) area the spectral density becomes non-analytic. However, the winding 
Wilson loops are perfectly analytical as emphasized by Narayanan and Neuberger [1], so 
the existence of the critical loop size does not represent a phase transition in these loops. 
As a matter of fact, a superficial inspection of the Wilson loops defined on winding curves 
does not show any sign of a critical area. So the question is what is the meaning of the 
non-analyticity in the spectral density? 

In this note we shall show that the critical area is indeed encoded in the Wilson loops 
with more than one winding. This encoding shows up as a smooth transition from different 
types of physics. At short distances the winding loops oscillate, and have a "non-stringy" 
behavior. For larger loops the string-like behavior is regained. We compare this behavior 
successfully to tunneling in ordinary non-relatistic quantum mechanics. It turns out that 
the critical size corresponds to the turning point, and the smaller sizes are analogous to the 
"classically allowed" region, whereas the confining larger loops correspond to the "classically 
forbidden" region. In view of the results obtained in ref. [1] a similar phenomenon most 
likely happens in four dimensions. 

We emphasize that the idea of tunneling from an intermediate state with no confinement 
to a confining state is certainly not new. Neuberger [4] has discussed this several years ago, 
and it is also mentioned in ref. [1]. Our new point is that in contrast to what one would 
naively think by inspection of the winding Wilson loops, these have encoded information on 
a tunneling transition with a turning point exactly at the critical area. 

In Section 2 we summarize the results on two-dimensional N = oo QCD. In Section 3 we 
develop asymptotic results for Wilson loops with a large number of windings, and in Section 
4 we connect these loops with the gap in the spectral density. In Section 5 we develop the 
tunneling picture in the WKB approximation, and make some conclusions. 

2 The gap and non-analyticity in the spectral density 

We start by giving an overview of the results previously obtained [2, 3] in two dimensional 
QCD which are relevant for the following. The spectral density gives the distribution of the 
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continuous N = oo eigenvalues 9. From [2], eqs. (4.1) and (4.2), one easily finds 
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where y = —tyApa(9), and where A = g 2 Nx(Area), which means that A = 2k, where k is 
the parameter used in [2]. The parameter A is a dimensionless area, defined such that the 
Wilson loop has the behavior 

W = e- A '\ 



There exists a finite gap 9 C < n for A < 4, given by 



/ A 2 A 
6 C = W A - — + arccos(l - — ). 

For A > 4 the gap disappears, and the maximum value of 9 is ir. 
Near the gap one has 
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Thus the derivative dpA{0)/dd becomes infinite at the gap. 
Close to A = 4 this breaks down, and instead one gets 
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For A> 4 the behavior near 9 = ir is smooth. The value y c at # = 7r is determined by 
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Expanding around y = y c with y m y c + 5y, we get 
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From this we see that the derivative dp/ 89 vanishes at 9 = it. For A large, y c —^4(1 — 
2 e _A / 2 )/2, and one obtains 



PA(^PA(7r) + -^e- A / 2 (7r-#) 2 + 
For A — > 4+ we can easily find the small y c from (6), 
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Therefore the derivative &pa{9) / dA becomes infinite for A — > 4. In the following we shall 
show that this behavior reflects itself in the winding Wilson loops as tunneling with the 
turning point A = 4. 
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3 The Wilson loops with windings 



The Wilson loop with n windings found by Kazakov and Kostov [7] is completely analytic 
in A, 

W^ n \A) = i Ll_ x {nA) e~ nA/2 , (10) 
where L\_ x is the Laguerre polynomial of type one, 
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The lack of analyticity of p for A = 4 obviously occurs due to the summation defining the 
spectral density, 

i r oo 

1 + 2]Tcos(n#) W {n \A) . (12) 

?1=1 



Also, from (10) it appears that has no reference whatsoever to the occurence of the 

gap 8 C in the spectral density. 

The non-analyticity and the existence of a gap for A < 4 is clearly related to the large n 
behavior of W^ n \ since if the sum in (12) is truncated, these phenomena do not occur. This 
is well known to several authors (for a recent reference, see [1]). However, we shall now show 
that W^ n \A) has a behavior which, at least for large n, reflects the existence of the gap 9 C . 
To this end we need the uniform large n asymptotic behavior of the Laguerre polynomial 
found by Erdelyi [5] from Laguerre's differential equation, and later also found from inte- 
gral representations of the Laguerre polynomial [6]. There are two (partially overlapping) 
representations of these asymptotes, called the Bessel form and the Airy form. Applying 
Erdelyi's Bessel form to eq. (10) we get 

W^\A) « ^ ( ^EZ±^^j 3l (2n(V^¥ + arcsinv^)) + ... , (13) 



where t = Aj '4. We have checked numerically that the corrections are extremely small 
even for moderate values of n, except for A w 4, where the approximation breaks down. 
However, the larger n is, the closer to A = 4 does the approximation work. 
The Bessel form is thus valid for A < 4 and A ^ 0. 

On the other hand, applying Erdelyi's Airy form, we obtain for < t < 1 



W {n >(A) ss '— — At ( - (3n) 3 (arccos Vt - Vt - t 2 ) s ) 



For t > 1 we get from [5] 
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The Airy function is defined by 

I c°° L 

Ai(z) = - / cos(-t 3 + zt)dt (16) 
n Jo 3 

and the Airy forms (14) and (15) can be expressed in terms of Bessel functions by means of 
the formulas 

Ai(z) = -y/\z\ f J 1/3 (-Li_ )+ J -i/3(^— J J forz<0, (17) 

and 

Ai(z) = tf 1/3 for * > 0. (18) 

It should be noticed that these Bessel functions are different from the J\ Bessel function in 
eq. (13). 

4 Connection between the gap parameter 6 C and 

In this Section we shall show that in spite of their appearence the W^^s with n > 1 know 
about the gap, and have quite different behavors on the two sides of A = 4. 

If we consider A < 4 we see that the argument of the Bessel function (13) is similar to 
the gap parameter 9 C given in eq. (3). The connection is easily found, 



\ft — t 2 + arcsin \/t = ^ A — + 2 arccos ^1 — ^ 
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Therefore eq. (13) can be written 

w '"'<*> = ^(4M-i)v« *W + -"- < 20 > 

Similarly, we also have 

arccos - Vt^¥ = hn - 6 C ), (21) 

so we can rewrite eq. (14), 

^"'^ k U^-i^ ^(-(|) 2/3 (" - <y 2/3 >- < 22 > 

By means of eq. (17) this can also be written in terms of Bessel functions, 

W^(A) « 1 J ; [Ji/ 3 (n(7r - C )) + J- 1/3 (n(ir - 6 C ))] . (23) 



The various Bessel functions in these expressions oscillate, and the oscillations are governed 
entirely by the gap parameter, either 9 C directly as in (20), or else by the length of the gap 
7r — 9 C as in (23). Thus, the Wilson loops with windings n > 1 actually know about the gap 
parameter. 

We have checked numerically that (20) and (23) give a very precise representation of the 
exact expression for W {n) in eq. (10). This holds for surprisingly low values of n, especially 
as far as the Airy form (23) is concerned. With the Bessel form (20) there are deviations 
close to A = 4, but one can get closer by taking large n's, e.g. n = 10. 

When n is sufficiently large, and 9 C or n — 9 C are not too small, we can use the well known 
asymptotic representations for the Bessel functions to rewrite (20) and (23), 
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Although the Bessel functions are different in eqs. (20) and (23), after use of the addition 
theorem for cosines they lead to the same expression (24), as one would expect. The asymp- 
totic expression (24) can be derived in a straightforward manner by the saddle point method 
from the integral representation of the Laguerre polynomial, 



-nAt 
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(25) 



where the contour C encloses the origin. 

The asymptotic expression (24) gives good results if A is not too close to 4. As an 
example, let us check the zeros predicted by (24) against the zeros in L^_ 1 (nA) for the 
rather moderate value n — 5. The 9 c (A)'s which correspond to the zeros in L\(5A) are 



0.766, 1.403, 2.034, and 2.662 
The zeros predicted by (24) are 

0.785, 1.414, 2.042, and 2.670. 



(26) 



(27) 



The largest discrepancy is of order three per cent. With n > 5 the accuracy improves, of 
course. 

We now turn to A > 4, where we have eq. (15). By use of (18) we can rewrite this 
expression as 



W {n \A) 
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By means of the asymptotic expansion of the if— function this gives (A > 4) 
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This result can also be derived by the saddle point method from the integral representation 
(25). For A -> oo eq. (29) leads to 



W^ n \A) w ^tl e~ nA /\ (30) 

V27T n 3 ' 2 



This is precisely the expression one would get from the exact formula (10) by taking into 
account only the first term in the Laguerre polynomial (11) and using Stirling's formula to 
estimate the factorials for n — > oo. 

It is easy to check that eqs.(23) and (29) give good numerical results for A > 4. The 
same is not true for (30), which is only valid for A's somewhat larger than 4. 

The turning point region around A = 4 can be approximated by simpler expressions. 
The results depend on how close A is to 4. If A — 4 << n~ 2 / 3 one finds [5] 

(_l)n-l 

w ' n> ^ M { 2„)'m/ 3 y (31 > 

whereas for A — 4 << n~ 2 / 5 the approximation is 

W™(A) « M{n^{A - 4)/(22 1 /3)). ( 32 ) 

These two expressions agree if A = 4 because Ai(0) = l/(3 2//3 r(2/3)). 



5 On the physical meaning of the gap and non-analyticity: 
Tunneling 

The results obtained in the previous section show that the Wilson loop with two or 

more windings is influenced by the gap and the A = 4 behavior. There are two different 
physical regions, but they are smoothly connected: In the region A < 4 there is no area 
behavior, but oscillations. In the second case, there is an area behavior. This is in contrast 
to W^-\ where there is always an area behavior e~ A l 2 . 

One may wonder why the area behavior is not present for A less than 4, because the 
Wilson loop (10) contains the factor e~ nA / 2 , which has a large negative exponent for large 
n, even if A is moderate. The answer to this point is that the Laguerre polynomial has large 
coefficients for large n. For example, the leading term is of order A n e n , which for moderate 
A can win over e~ nA l 2 . 

The non-analytic behavior of the spectral density must be a non-smooth reflection of the 
smooth transition between the two physically different regions, separated by A = 4. The 
Wilson loops are analytic in A, as follows trivially from (10). However, if we think in 
terms of a string description of the dynamics, we see that the Wilson loop with n > 1 is able 
to go from the "evident" string at larger A to something which does not look "stringy" at 
smaller A's. 

Let us look at the eigenvalue gap parameter C (A) as a coordinate, which is compact but 
continuous (thanks to iV = oo) in the region A < 4. For A > 4 we can instead take the 
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argument of the Bessel i^-function as a coordinate. Let us denote this set of "coordinates" 
by & 



£i=tt-0 c (A) (^<4), 6 = V A V4 - A - 2 ln( v^/2 + ^4/4 - 1) (A > 4). (33) 

Then the (smooth) transition that occurs at A = 4 looks very much like tunneling, where 
A < 4 is the "classically allowed region", whereas A > 4 is the region which is "classically 
forbidden", and A = 4 is the turning point. The Bessel functions we found are indeed similar 
to those wave functions that are encountered in the WKB approximation. 

To see the connection to the WKB approximation in details, we notice that £ should be 
represented by an integral over "momenta" ki and kit" . The integrals are 

rA rA 

/ k^dx — n £i and / k R dx = n £ 2? (34) 

J A J 4 

where we consider n to be a fixed number. From eqs. (33) we get 




n A — A /0 . 

k * = 2 v— < 35) 

Near the turning point A = 4 we therefore see that the squares of the momenta are linear, 
oc 4 — A or A — 4. Therefore the WKB-solution on the left should be expresed in terms of 
$ 2 J±i/z{n^i) and to the right by £2^^1/3(^2) (see e.g. ref. [8], pp. 188-190). Taking into 
account that 

%(x) = ^[L 1/3 (x)-/ 1/3 (x)]. (36) 

this is precisely what we obtained using Erdelyi's approximations. This even applies to the 
prefactors 

(4/A-l)- 1 / 4 or (1 -4/A)- 1 / 4 , (37) 

which correspond to the prefactors k^ 2 or k^ 2 occuring in the WKB expressions [8]. 

It should be noticed that even without using the large n approximate results, the Wilson 
loops (10) for n > 1 clearly distiguish between values of A larger or smaller than 4. This for 
example follows because the Laguerre polynomials L°_ 1 (x) have all their zeros in the region 
< x < v, where v = 4n + 2a — 2, and outside this interval the polynomials are monotonic 
[5]. In our case (take x = nA) this condition for oscillations means < A < 4, so the zeros 
are all below A = 4, and in this region the Wilson loops (10) oscillate. There is therefore 
always a non-stringy behavior below A = 4 for n > 2. In this connection it is amusing to 
notice that if we define a "wave function" ip 

ij n (A) = AW (n \A), n>l, (38) 

then from (10) and the differential equation for the Laguerre polynomials it follows that ip 
satisfies the linear "Schrodinger equation" 



dA 2 4 \A 
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In this equation no approximation has been used. Eq. (39) clearly shows the existence of 
the turning point at A = 4. The approximate results discussed in the preceeding sections 
were actually derived from (39), see ref. [5]. 

The physics behind the domains separated by A = 4 is indicated by the spectral density. 
In the extreme case where A is small, the eigenvalues 6 are centered around 9 = 0, and in 
the other extreme where A is large, the eigenvalues are almost uniformly distributed. Thus, 
the domain to the right of A = 4 corresponds approximately to a completely disordered 
state, whereas the domain to the left of A = 4 is a much more ordered state, at least as far 
as the eigenvalues are concerned. Therefore the confinement is related to randomness, and 
the approximately ordered state leads to an "evaporation" of the effective string tension. 
This evaporation happens to the left and even somewhat to the right of of A = 4, see the 
exponent in eq. (29). 

The coordinates £1,2 represent a new degree of freedom, which perhaps may be interpreted 
as an extra dimension, following ideas by Das and Jevicki [9] and by Brower and collaborators 
[10]. What this adds up to is that £i )2 might play a dynamical role in an effective string 
theory in three dimensions, where the action provides a barrier, which produces tunneling 
in the new dimension. One difficulty with this picture is that this tunneling does not occur 
for W^. 

6 Conclusions 

The main conclusion of this note is that the behavior of W^ n \A) as a function of A physi- 
cally represents a smooth tunneling between a short distance (A < 4) and a large distance 
state (A > 4). In view of the results obtained by Narayanan and Neuberger [1] on the 
relation between four and two dimensional QCD this may be true not only in two, but 
also in four dimensions. The A < 4 region may be considered as a transient state where 
there is no "obvious" string-picture. The tunneling to the disordered confined state has un- 
doubtedly different causes in different dimensions (instantons, monopoles, vortex lines,...). 
In four dimensions the asymptotic freedom short-distance region would be hooked on to the 
approximately ordered state, which is then intermediate between freedom and confinement. 
Of course, such a picture does not exclude that there exists a sufficiently intelligent string 
theory which can accomodate this behavior, as advocated by several authors, most recently 
by Brower [10]. However, it is of some interest that the W^ n ^s have at least some knowledge 
of the complexity involved at different scales. 

The connection between disorder and large distance dimensional reduction was actually 
argued a long time ago [11] to be connected to confinement, because e.g. four dimensional 
QCD could reduce to two dimensional QCD at large distances as a consequence of the 
disorder. The results obtained in [1] precisely means the existence of a disordered vacuum 
at large distances in four dimensions. The phenomenon disorder — > reduced dimension 
is known to occur in simple condensed matter models [12], e.g. a scalar field coupled to 
a random Gaussian field. In this relatively simple case the disorder reduces the effective 
dimension D by two, D — > D — 2. The complexity of QCD means that the results and the 
methods from solid state physics have not been generalized to QCD, and it is not known 
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by how much the dimension effectively reduces. Thus it may be that both four and three 
dimensions effectively reduce to two dimensions. In QCD the disordered field should be 
generated dynamically in the "true" vacuum, and not put in externally, as in the condensed 
matter case. The complexity of QCD is a good reason to study this problem by using lattice 
QCD. It would be interesting to investigate on a lattice whether the tunneling transition 
found in the present paper also occurs for (winding) Wilson loops in four dimensions. The 
similarity of the spectral density in four and two dimensions found in [1] indicates that this 
would be the case. 
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